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1.  EXECUTIVE  SUMMARY  * 


1.1.  Overview  of  Project 

While  it  is  generally  accepted  that  near  source  nonlinearities  in  rock  deformation 
influence  the  far-field  seismic  signature  of  a  seismic  source,  an  assessment  of  the 
importance  of  these  nonlinearities  has  been  hampered  by  a  lack  of  well  constrained 
rheologies  with  which  to  perform  forward  modeling  simulations.  Previous  work  aimed 
at  developing  nonlinear,  hysteretic  rheologies  have  largely  been  based  on  experiments 
covering  limited  stress  conditions  and  modes  of  deformation.  While  these  experiments 
have  been  important  in  guiding  model  development,  the  resulting  models  have  been 
limited  in  their  value  due  to  the  restricted  loading  paths  for  which  they  are  valid. 

In  order  to  develop  forward  models  for  wave  propagation  in  the  moderate  strain 
and  damage  regimes,  experiments  using  a  much  wider  range  of  loading  path  were 
needed  to  test  and  constrain  models  suitable  for  arbitrary  loading  paths.  In  this  study, 
an  experimental  protocol  was  developed  and  implemented  which  provides  detailed  con¬ 
straint  of  the  nonlinear,  hysteretic  deformation  of  rock  subject  to  arbitrary  perturbations 
in  axially  symmetric  compressive  stresses.  Experiments  were  performed  on  a  diverse 
suite  of  samples,  including  Berea  sandstone,  Sierra  White  granite,  and  tuff  from  N- 
tunnel  at  the  Nevada  Test  Site  (NTS). 

1.2.  Project  Goals 

Prior  to  the  start  of  this  project,  a  rheological  model  had  been  developed  for  the 
case  of  uniaxial  stress  (see  Boitnott,  1993).  This  model  was  thought  to  be  illustrative 
of  the  fundamental  characteristics  of  nonlinearities  in  rock  deformation  in  the  moderate 
strain  regime,  but  was  limited  in  its  applicability  due  to  the  restrictive  loading  path  for 
which  it  was  valid.  The  underlying  goal  of  this  project  was  to  extend  this  rheology  to 
be  applicable  for  more  general  loading  paths. 

In  order  to  apply  the  developed  rheology  directly  to  problems  of  wave  propaga¬ 
tion  from  a  seismic  source,  the  rheology  needed  to  be  extended  to  loading  paths 
expected  in  the  field.  To  accomplish  this  task,  four  areas  needed  to  be  addressed. 
First,  the  rheology  had  to  be  extended  to  include  a  wide  variety  of  loading  conditions. 
In  order  to  model  compressional  wave  propagation  (including  the  geometric  effects  of 
spherical  wavefronts)  the  uniaxial-stress  rheology  needed  to  be  extended  to  loading 
paths  which  can  be  transitional  between  uniaxial-strain  and  uniaxial-stress.  In  addition, 
in  order  to  model  shear  wave  propagation,  laboratory  data  was  needed  to  constrain 
hysteresis  in  the  shear  modulus  during  shear  loading.  Second,  the  effects  of  confining 
pressure  on  the  rheological  models  needed  to  be  constrained  in  order  to  address  ques¬ 
tions  concerning  the  effect  of  depth  of  burial  on  the  seismic  signature  of  an  under¬ 
ground  explosion.  Third,  data  was  needed  on  tuff  in  order  to  make  use  of  well  cali¬ 
brated  field  data  from  the  Nevada  Test  Site.  Fourth,  the  rheology  needed  to  be 
extended  into  the  brittle-damage  regime. 

1.3  Summary  of  Findings 

Uniaxial  stress  experiments  similar  to  those  performed  previously  (but  in  this 
study  performed  as  a  function  of  confining  pressure)  confirmed  that  the  uniaxial  stress 
rheology  simplifies  with  increasing  confining  pressure.  Deformation  becomes  more 
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linear  and  less  hysteretic  with  increasing  confining  pressure,  although  the  observed 
nonlinearities  clearly  persist  at  all  stress  conditions  studied.  Complex  rheological  pro¬ 
perties  at  reversals  in  loading  direction  (i.e.  A E  effect  ,  see  Boitnott,  [1993])  were 
found  to  be  less  evident  at  elevated  confining  pressures.  As  a  result,  the  eventual 
rheology  was  simplified  to  make  use  of  a  common  parameterization,  denoted  the 
"asymmetric  bow-tie  parameterization"  (ABP),  which  appears  well  suited  for  modeling 
hysteresis  in  stress  versus  strain  for  a  wide  variety  of  loading  paths.  It  is  found  that 
the  ABP  produces  reasonable  fits  to  data  for  all  loading  paths  and  at  all  stress  condi¬ 
tions  studied,  provided  that  the  stress  perturbations  are  not  exceptionally  large  and  the 
loading  perturbation  is  restricted  to  a  single  mode  of  deformation  (see  section  3). 

Radial  strain  measured  during  the  uniaxial  stress  experiments  provided  a  test  of 
whether  a  simple  rheology  could  be  constructed  by  assuming  that  the  dominant  source 
of  hysteresis  was  in  the  shear  component  of  deformation.  This  was  found  not  to  be  the 
case,  as  the  experiments  indicated  appreciable  hysteresis  (and  thus  energy  loss)  in  the 
volumetric  strain  versus  mean  stress  relationship  under  some  conditions.  Preliminary 
observations  of  radial  strain  during  uniaxial  stress  perturbations  lead  to  an  initial 
attempt  to  modify  of  the  uniaxial  stress  rheology  to  include  the  effects  of  changing 
confining  pressure  (section  2.2).  While  this  attempt  produced  realistic  stress-strain 
curves  for  uniaxial  strain  paths,  the  model  was  clearly  an  over-simplification  as  it  did 
not  hold  up  well  when  used  to  predict  uniaxial  stress  data  as  a  function  of  confining 
pressure  (section  2.3).  As  a  result,  development  of  a  rheology  based  on  uniaxial  stress 
and  confining  pressure  was  abandoned  in  favor  of  a  different  rheological  description 
based  on  the  decomposition  of  stresses  into  mean  stress  and  shear  stress  components. 

Hydrostatic  stress  cycling  experiments  confirmed  that  the  hysteresis  in  the  mean 
stress  versus  volumetric  strain  relationship  was  appreciable,  and  that  it  could  be 
characterized  by  a  qualitatively  similar  rheology  as  previously  developed  for  the  uniax¬ 
ial  stress  versus  axial  strain  data  (i.e.  can  be  described  with  an  ABP  for  limited  stress 
perturbations).  Similarly,  torsion  experiments  on  solid  cylinders  were  used  to  explore 
the  rheology  for  the  shear  deformation.  Again  the  rheology  was  found  to  be  describ- 
able  using  the  ABP. 

The  observations  that  the  mean  stress,  shear  stress,  and  uniaxial  stress  data  could 
all  be  modeled  reasonably  well  with  the  same  parameterization  lead  to  the  idea  that  the 
uniaxial  stress  data  might  be  modelable  using  a  superposition  of  the  mean  stress  and 
shear  stress  rheologies.  However,  it  was  found  that  a  simple  superposition  was  not  sup¬ 
ported  by  the  data,  as  the  parameters  needed  to  fit  the  mean  stress  versus  volumetric 
strain  relationship  during  uniaxial  stress  perturbations  were  not  quantitatively  consistent 
with  those  predicted  from  pure  mean  stress  experiments.  Thus,  the  superposition  was 
not  directly  generalizable  to  arbitrary  loading  paths,  indicating  that  "coupling"  had  to 
exist  between  the  mean  stress  and  shear  stress  rheologies. 

The  coupling  terms  were  clearly  complex  enough  that  it  was  not  viable  to 
develop  a  rheology  based  on  combining  results  from  separate  experiments  characteriz¬ 
ing  individual  parts  of  the  rheology.  This  lead  to  the  development  of  a  testing  protocol 
which  would  have  sufficient  flexibility  to  constrain  a  more  general  rheology,  including 
mixed  mode  loading  paths.  An  experimental  protocol  was  developed  and  implemented 
to  constrain  the  nonlinear,  hysteretic  deformation  of  rock  in  the  moderate  strain 
regime.  While  more  work  is  required  to  fully  integrate  the  coupling  terms  into  a 
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hysteretic  rheology,  the  resulting  rheology  has  been  developed  to  a  point  where  for¬ 
ward  modeling  of  compressional  wave  propagation  can  be  performed  in  order  to  assess 
which  features  in  the  rheology  are  important  to  a  particular  application. 

Tests  using  loading  paths  which  quasi-statically  mimic  the  loading  path  expected 
from  a  spherically  propagating  compressional  wave  were  performed  to  assure  that  the 
developed  rheology  was  capable  of  handling  explosion  like  loading  paths.  It  is  found 
that  for  a  broad  class  of  explosion  like  loading  paths,  a  simpler  rheology,  using  an 
ABP  constrained  from  uniaxial  strain  data,  may  be  sufficient  for  many  applications. 

In  order  to  test  the  rheology  further  and  explore  the  physical  meaning  of  some  of 
the  parameters,  tests  were  designed  and  executed  which  combined  the  static  stress 
versus  strain  response  with  measurements  of  ultrasonic  velocities.  The  velocities  were 
used  to  infer  the  small  strain  dynamic  moduli,  and  those  values  were  compared  with 
the  moderate  strain  rheology.  The  results  indicated  that  the  moderate  strain  rheology 
does  have  the  proper  behavior  in  the  limit  of  small  strains,  as  it  approaches  a  linear 
elastic  model  consistent  with  the  measured  velocities. 

Exploratory  tests  were  also  performed  on  Barre  granite  to  probe  the  transition  into 
the  damage  regime.  The  experiments  illustrated  the  added  complexity  of  experimental 
characterization  in  the  damage  regime  due  to  time  dependent  deformation.  They  also 
provided  some  indication  that  pressure  effects  on  the  uniaxial  stress  rheology  are  simi¬ 
lar  both  within  the  moderate  strain  and  damage  regimes  (section  5). 

1.4  Collaborative  Efforts 

One  emphasis  of  this  project  was  to  support  other  work  within  the  AFOSR  and 
related  programs  by  providing  laboratory  data  to  others  for  the  constraint  of  modeling 
efforts  of  near  source  processes  and  wave-propagation.  Experiments  performed  during 
the  course  of  this  project  were  continually  revised  in  order  to  better  serve  these  pro¬ 
jects.  Projects  specifically  targeted  for  support  involved:  the  development  of  an  endo- 
chronic  rheology  for  modeling  wave  propagation  in  nonlinear  materials  by  J.  B.  Min¬ 
ster  (U.  Cal.  San  Diego)  and  S.  Day  (San  Diego  State  Univ.);  development  of  a  dam¬ 
age  mechanics  model  for  wave  propagation  in  the  damage  regime  by  C.  Sammis 
(USC);  and  development  of  an  "equivalent  elastic"  method  for  modeling  wave  propa¬ 
gation  through  nonlinear  materials  by  L.  Johnson  (U.  Cal.  Berkeley  ).  On  two  occa¬ 
sions,  students  working  on  these  projects  visited  NER  and  participated  in  the  design, 
execution,  and  preliminary  analysis  of  laboratory  experiments. 

Inna  Altschul  came  from  USC  to  work  on  tests  to  help  characterize  rock  rheology 
in  the  damage  regime.  A  number  of  tests  were  performed,  and  the  data  became  a  cen¬ 
tral  part  of  her  masters  thesis  [Altschul,  1995].  H.  Xu  came  from  U.  Cal.  San  Diego 
to  work  on  extending  the  endochronic  rheology  previously  established  by  Day  et  al. 
[1993].  Specific  experiments  were  designed  to  constrain  the  hysteretic  rheology  during 
arbitrary  loading  paths  involving  simultaneous  changes  in  both  axial  and  confining 
stresses.  The  results  have  subsequently  been  used  for  wave  propagation  modeling  [Xu 
et  al.  1996].  The  data  and  resulting  model  will  continue  to  be  used  as  part  of  Xu’s 
PhD.  thesis.  Experiments  were  also  performed  specifically  to  support  Lane  Johnson’s 
development  of  the  equivalent  elastic  method  for  modeling  wave  propagation  (see 
Johnson  [1996]),  with  the  idea  that  laboratory  data  on  tuff  was  needed  to  help  con¬ 
strain  interpretation  of  field  data  from  the  NPE  test  at  NTS.  While  the  data  has  not  yet 
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been  used  directly  in  this  context,  data  on  tuff  from  N-tunnel  at  NTS  is  now  available 
to  help  constrain  forward  models  of  NTS  explosions. 

Other  researchers  also  made  use  of  the  data  set  for  the  purpose  of  constraining 
models  for  nonlinear  deformation  of  rocks.  D.  Lieblich,  a  PhD.  student  from  U. 
Conn.,  was  provided  data  to  support  the  development  of  an  adhesion  based 
micromechanical  model  for  hysteresis  in  rock  deformation.  Data  was  also  provided  to 
K.  McCall  (U.  Nev.  Reno)  and  R.  Guyer  (U.  Mass.  Amherst)  for  use  in  further 
development  of  a  discrete  element  model  for  nonlinear  deformation  of  rocks. 


1.5  List  of  Publications 

The  following  publications  resulted,  either  in  whole  or  in  part,  from  this  project.  Com¬ 
plete  citations  can  be  found  in  section  6. 

Altschul,  I.,  The  relationship  of  stress  to  strain  in  the  damage  regime  for  a  brittle  solid 
under  compression,  M.S.  Thesis,  University  of  Southern  California,  71  pp.,  1995. 

Boitnott,  G.  N.,  Experimental  characterization  of  the  nonlinear  rheology  of  rock,  PL- 
TR-95-2108,  1995. 

Boitnott,  G.  N.,  Constructing  a  general  rheological  model  for  rock  deformation,  PL- 
TR-96-2153,  p.  149-158,  1996. 

Boitnott,  G.  N.,  Experimental  characterization  of  the  nonlinear  rheology  of  rock,  Int. 
J.  Rock  Mech.  &  Min.  Sci.,  34:3-4,  Paper  No.  033,  1997. 

Guyer,  R.  A.,  K.  R.  McCall,  G.  N.  Boitnott,  L.  B.  Hilbert  Jr.,  and  T.  J.  Plona,  Quanti¬ 
tative  implementation  of  Preisach-Mayergoyz  space  to  find  static  and  dynamic 
elastic  moduli  in  rock,  J.  Geophys.  Res.,  v.  102,  n.  B3,  p.  5281-5293,  1997. 

Sammis,  C.,  Stress  strain  modeling  in  the  damage  regime,  PL-TR-95-2108,  p.  491-495, 
1995. 

Xu,  H.,  S.  M.  Day,  and  J-B.  H.  Minster,  Two-Dimensional  nonlinear  wave  propagation 
in  an  endochronic  material,  PL-TR-96-2153,  p.  301-308,  1996. 
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2.  EXPLORATORY  EXPERIMENTS 


2.1  Overview 

The  first  task  was  to  design  experiments  which  would  be  suitable  for  characteriz¬ 
ing  the  rheology.  As  the  nature  of  the  rheology  was  largely  unknown,  a  number  of 
exploratory  experiments  were  performed  to  develop  baseline  observations  from  which 
to  identify  the  critical  elements  required  for  a  meaningful  rheology.  In  this  section  we 
review  these  experiments  and  provide  examples  of  results  of  each  type.  The  results  are 
discussed  in  the  context  of  how  they  contributed  to  the  development  of  the  rheological 
model  described  in  section  3. 

2.2  Initial  Attempts  to  Extend  the  Uniaxial  Stress  Rheology 

As  a  starting  point  to  extend  the  rheology  to  loading  conditions  other  than  uniax¬ 
ial  stress,  unconfined  uniaxial  stress  experiments  were  performed  on  Sierra  White  gran¬ 
ite  during  which  both  axial  and  radial  strain  were  measured.  These  experiments  were 
motivated  by  the  need  to  identify  simplifying  principles  which  might  be  used  to  extend 
the  uniaxial-stress  rheology  to  the  case  of  arbitrary  (triaxial)  stress  perturbations. 

An  example  of  data  from  an  experiment  on  Sierra  White  granite  is  shown  in  Fig¬ 
ure  1.  In  Figure  la  we  see  the  familiar  cusped  hysteresis  loop  between  axial  stress  and 
axial  strain  for  a  cyclic  axial  load  in  unconfined  compression.  In  Figure  lb,  we  plot  the 
corresponding  radial  strain  as  a  function  of  axial  stress.  Note  that  the  radial  strain  exhi¬ 
bits  little  measurable  hysteresis.  Importantly,  this  observation  rules  out  the  possibility 
of  decomposing  an  arbitrary  deformation  into  a  non-hysteretic  dilation  and  hysteretic 
distortion,  since  the  data  indicates  considerable  hysteresis  in  the  volumetric  strain  dur¬ 
ing  perturbations  in  uniaxial  stress. 

The  relative  lack  of  hysteresis  in  the  radial  strain  during  unconfined  compression 
did  however  indicate  a  possible  simplifying  principle  which  could  lead  to  a  more  gen¬ 
eral  rheology  capable  of  handling  loading  paths  transitional  between  uniaxial-strain  and 
uniaxial-stress.  The  data  in  Figure  lb  suggested  that  the  radial  (Poisson)  expansion 
was  a  dominantly  elastic  phenomena.  It  would  follow  that  the  confining  pressure 
required  to  maintain  uniaxial-strain  conditions  should  also  be  relatively  non-hysteretic 
with  respect  to  axial  load.  Based  on  this  observation,  we  initially  hoped  to  develop  a 
rheological  model  which  involved  the  addition  of  a  uni-valued  (non-hysteretic)  depen¬ 
dence  of  the  axial-strain  vs.  axial-stress  relationship  on  the  hydrostatic  component  of 
the  stress  state. 

The  notion  that  the  difference  between  uniaxial-strain  and  uniaxial-stress  loading 
may  be  included  in  the  rheology  with  the  rather  simple  addition  of  a  non-hysteretic 
confining  pressure  effect  was  consistent  with  some  limited  observations.  In  Figure  2 
we  compare  experimental  data  on  Berea  sandstone  from  Hilbert  et  al.  [1994]  with  a 
simple  extension  of  a  rheological  model  previously  developed  for  the  case  of 
unconfined  compression  (i.e.  from  Boitnott ,  [1993])  .  Figure  2a  shows  a  plot  of  axial 
strain  versus  axial  stress  for  uniaxial-stress  and  uniaxial-strain  loading  paths  on  a  sam¬ 
ple  of  Berea  sandstone.  Model  predictions  based  on  this  extension  of  the  uniaxial 
stress  rheology  are  shown  in  Figure  2b.  The  model  results  for  uniaxial-strain  have 
been  produced  by  adding  a  confining  pressure  effect  to  the  uniaxial-stress  model.  We 
start  by  defining  the  axial  stiffness  r|  =  don/dej],  where  Gn  and  en  are  the  axial 
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stress  and  axial  strain  respectively.  We  then  assume  that  r \  =  XE,  where  E  is  the  hys- 
teretic  Youngs  modulus  describing  the  rheology  for  uniaxial-stress  perturbations,  and  X 
is  a  scaling  factor  which  is  a  linear  function  of  the  confining  stress.  The  confining 
stress  is  assumed  to  change  linearly  with  axial  stress  based  on  the  observed  lack  of 
hysteresis  in  the  radial-strain  versus  axial-stress  relationship  as  discussed  in  the  context 
of  Figure  lb.  Note  that  the  model  results  compare  favorably  with  the  experimental 
data,  illustrating  that  the  scaled  rheology  is  reasonable. 

2.3  Explicit  Tests  of  the  Effects  of  Confining  Pressure 

Explicit  tests  of  this  scaled  rheology  were  made  by  performing  uniaxial  stress 
tests  similar  to  those  performed  previously,  but  with  the  addition  of  confining  pressure 
as  a  test  variable.  Tests  were  performed  on  Berea  sandstone,  Sierra  White  granite,  and 
tuff  from  NTS. 

An  example  of  the  loading  protocol  used  on  Berea  Sandstone  is  shown  in  Figure 

3.  Stress  versus  strain  response  for  two  of  the  loading  sequences  are  shown  in  Figure 

4.  Axial  stress  versus  axial  strain  response  at  pressure  is  qualitatively  similar  to  that 
observed  for  the  unconfined  tests.  The  same  parameterization  as  used  in  Boitnott 
[1993]  is  found  to  be  suitable  at  elevated  confining  pressure.  As  expected,  the  rheol¬ 
ogy  does  appear  to  simplify  somewhat  with  pressure,  with  the  A E  term  getting  smaller 
(particularly  noticeable  for  the  transition  from  loading  to  unloading)  and  the  hysteresis 
loops  getting  stiffer,  straighter,  and  more  symmetric. 

Attempts  to  fit  the  pressure  effect  with  a  simple  X  term  as  described  above  were 
not  very  successful,  particularly  for  Berea  sandstone.  Some  representative  plots  of 
Young’s  modulus  versus  differential  stress  at  various  confining  pressures  are  shown  in 
Figures  5-7.  It  is  clear  from  the  data  that  in  general,  X  must  be  a  function  of  both 
confining  pressure  and  axial  load.  The  effect  of  confining  pressure  (i.e.  effective  X.) 
appears  larger  at  low  axial  stress,  and  decreases  with  increasing  axial  stress.  More 
problematic  is  the  fact  that  X  for  loading  does  not  appear  consistent  with  the  required 
X  for  unloading  (even  for  a  fairly  restricted  stress  range).  These  discrepancies  are  par¬ 
ticularly  true  for  Berea,  but  also  present  to  some  extent  in  both  the  granite  and  tuff 
data.  Thus  the  notion  of  X  as  a  rheological  parameter  was  abandoned,  although  it  still 
might  be  considered  suitable  in  some  cases  as  an  approximate  rheology  (see  for  exam¬ 
ple  section  5). 

Other  observations  lead  to  the  abandonment  of  the  uniaxial  stress  and  confining 
stress  parameterization  altogether.  Unlike  the  data  for  granite  discussed  in  section  2.2, 
significant  hysteresis  is  observed  between  the  radial  strain  and  the  axial  stress  for 
Berea  and  tuff,  undermining  the  original  justification  for  X  (see  Figures  8b,  9b,  10b). 
Hysteresis  between  the  radial  strain  and  axial  stress  is  also  somewhat  apparent  in  the 
granite  at  elevated  confining  pressures,  indicating  a  pressure  dependence  of  this  hys¬ 
teresis. 

The  failure  of  the  uniaxial  stress  based  model  lead  to  the  exploration  of  a  rheol¬ 
ogy  based  on  the  mean  stress  and  the  shear  stress  (rather  than  the  confining  pressure 
and  the  differential  stress).  Data  from  experiments  such  as  that  in  Figure  3  were 
decomposed  into  mean  stress,  shear  stress,  volumetric  strain,  and  shear  strain.  This 
view  of  the  data  is  plotted  in  Figures  8-10(c,d).  This  data  seemed  to  indicate  that  a 
rheology  parameterized  in  this  way  might  be  more  flexible,  and  could  still  build  on  a 
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parameterization  like  that  used  previously  for  uniaxial  stress.  Both  the  mean  stress  and 
the  shear  stress  rheologies  appear  rather  simple,  with  hysteresis  loop  shapes  mimicing 
the  previously  observed  uniaxial  stress  behavior  in  many  ways  (note  the  similarity  in 
the  stress  strain  relationships  in  Figures  8-10:a,c,d).  Unlike  the  earlier  granite  data,  the 
data  on  Berea  (and  to  some  extent,  this  new  data  on  granite  as  well)  exhibited  little 
hysteresis  between  mean  stress  and  volumetric  strain  during  the  uniaxial  stress  cycles 
at  pressure,  while  tuff  exhibited  hysteresis  in  this  relationship  (see  Figure  8c,  9c,  10c). 
This  indicated  that  the  hysteresis  observed  during  uniaxial  stress  cycles  is  partitioned 
differently  for  different  rock  types.  This  partitioning  is  also  a  function  stress  condi¬ 
tions  and  loading  path  for  a  given  rock  type. 

2.4  Mean  Stress  (Hydrostatic)  Experiments 

The  observations  from  the  uniaxial  stress  tests  lead  to  the  need  for  some  baseline 
measurements  of  the  two  proposed  end-member  rheologies;  namely  those  for  pure 
mean  stress  perturbations  and  pure  shear  stress  perturbations,  respectively.  To  begin 
exploring  the  possibility  of  a  rheology  based  in  part  on  a  hysteretic  bulk  modulus,  a 
set  of  experiments  were  performed  to  characterize  the  rheology  describing  hydrostatic 
stress  and  volumetric  strain.  An  example  loading  protocol  is  shown  in  Figure  11. 

Each  of  the  samples  tested  exhibited  considerable  hysteresis  in  the  volumetric 
strain  versus  mean  stress  relationship.  Example  stress  strain  curves  are  shown  in  Fig¬ 
ures  12-14.  Comparing  Figures  12-14  with  the  mean  stress  plots  from  uniaxial  stress 
paths  in  Figures  8-10,  we  find  that  the  mean  stress  versus  volumetric  strain  relation¬ 
ships  are  different  for  uniaxial  stress  and  hydrostatic  compression.  For  each  rock-type, 
the  pure  hydrostatic  cycles  tend  to  exhibit  measurably  more  compliant  bulk  moduli 
than  observed  during  uniaxial  stress  cycles  for  similar  mean  stress  variations.  There 
also  appears  to  be  considerably  more  hysteresis  during  hydrostatic  cycling  than  during 
uniaxial  stress  cycling  (particularly  for  Berea  sandstone  and  Sierra  White  granite)  From 
these  observations  it  was  clear  that  observations  of  the  mean  stress  versus  volumetric 
strain  relationship  for  different  loading  paths  needed  to  be  made  in  a  more  integrated 
fashion  (i.e.  during  a  single  experiment)  in  order  to  constrain  any  systematics  in  the 
mean  stress  versus  volumetric  strain  relationship  for  different  loading  paths.  These 
experiments  are  discussed  in  section  3. 

2.5  Torsion  Experiments 

In  order  to  develop  similar  observations  for  shear  stress  behavior,  torsion  experi¬ 
ments  were  conducted  to  constrain  hysteresis  in  the  shear  modulus  as  a  function  of 
normal  load  and  shear  stress  history.  Experiments  were  performed  on  hollow  cylinders 
of  intact  Berea  sandstone.  Granite  and  tuff  were  not  tested,  as  concurrent  work  indi¬ 
cated  that  other  tests  were  more  easily  integrated  with  other  data. 

Examples  of  data  on  Berea  sandstone  are  shown  in  Figure  15.  In  Figure  15a,  the 
shear-stress  vs.  shear-strain  relationship  is  plotted  for  a  cyclic  torsional  load  at  a  fixed 
axial  load  of  5.4  MPa.  The  frequency  of  the  loading  oscillation  was  0.01  Hz  and  data 
from  two  different  amplitude  oscillations  are  shown.  Note  the  nested  nature  of  the 
hysteresis  loops  with  amplitude  and  the  fact  that  the  mean  slope  of  the  hysteresis  loops 
(i.e.  the  average  shear  modulus)  decreases  with  loop  amplitude.  Figure  15b  shows  the 
results  of  three  torsional  oscillations  at  three  different  axial  loads  (3.6,  5.4,  and  7.2 
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MPa).  As  expected,  the  average  shear  modulus  increases,  and  the  width  of  the  hys¬ 
teresis  loops  decreases,  with  increasing  axial  load. 

As  we  might  expect,  the  hysteresis  loops  appear  to  exhibit  the  same  "discrete 
memory"  and  "closed-loop"  characteristics  as  noted  in  studies  of  hysteresis  in  Youngs 
modulus  during  uniaxial-stress  perturbations.  Based  on  this  preliminary  data,  the  hys¬ 
teresis  in  the  shear  modulus  appears  relatively  easy  to  model  with  only  a  few  parame¬ 
ters.  As  an  example,  a  simple  "bow-tie"  model  for  the  shear  modulus  as  a  function  of 
shear  stress  represents  the  observations  fairly  well  (see  Figure  16).  In  Figure  17  we 
plot  the  predicted  stress-strain  behavior  for  the  experimental  cases  shown  in  Figures 
15a  and  15b. 

Comparing  these  results  to  those  from  the  uniaxial  stress  experiments,  general 
qualitative  agreement  was  found.  Both  the  shear  stiffening  with  normal  load  and  the 
decreased  width  of  the  hysteresis  loops  are  seen  in  the  torsion  and  uniaxial  stress  data. 
This  suggested  that  the  shear  rheology  may  be  easily  transportable  between  different 
loading  paths. 


3.  CONSTRUCTING  A  GENERAL  RHEOLOGICAL  MODEL 

3.1  Overview 

The  uniaxial  stress,  mean  stress,  and  torsion  experiments  were  exploratory  experi¬ 
ments  designed  to  develop  the  necessary  understanding  to  guide  development  of  a  real¬ 
istic  rheology.  Based  on  these  tests,  the  mean  stress  and  shear  stress  were  identified  as 
good  candidates  for  the  basis  of  a  general  rheology,  but  it  was  clear  that  coupling 
terms  which  link  shear  stress  effects  on  volumetric  strain  and  mean  stress  effects  on 
shear  strain  were  significant  enough  that  new  experiments  needed  to  be  designed  to 
support  such  a  rheology. 

Each  of  the  experiments  discussed  in  this  section  were  performed  on  all  three 
rock  types,  but  detailed  analysis  of  the  coupling  terms  was  restricted  to  data  on  Berea. 
Berea  sandstone  was  chosen  for  concentrated  effort  because  it  was  the  most  nonlinear 
of  the  three  rocks  being  studied.  In  this  section,  we  summarize  the  results  of  the 
experiments  and  discuss  how  they  are  used  to  develop  and  constrain  a  rheological 
model  for  nonlinear,  hysteretic  deformation  during  arbitrary  loading  protocols. 

3.2  Experimental  Protocol 

The  loading  history  of  an  exemplary  experiment  is  illustrated  in  Figure  18.  The 
experiment  involves  a  sequence  of  perturbations  in  load  subject  to  varying  boundary 
conditions  and  executed  from  a  series  of  initial  stress  conditions.  At  three  different 
starting  stress  conditions,  the  axial  stress  and  confining  pressure  are  controlled  to  pro¬ 
duce  perturbations  in  stresses  enforcing  the  following  boundary  conditions:  mean  stress 
perturbations  at  constant  differential  stress  (4  cycles  of  two  different  amplitudes), 
uniaxial  stress  perturbations  (2  cycles),  shear  stress  perturbations  at  constant  mean 
stress  (4  cycles  of  two  different  amplitudes),  and  uniaxial  strain  perturbations  (4  cycles 
of  two  different  amplitudes).  Detail  of  one  such  sequence  is  shown  in  Figure  19. 
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3.3  Formulating  a  Model 

The  goal  is  to  develop  a  model  capable  of  predicting  the  stress-strain  relationship 
for  arbitrary  loading  paths  in  the  stress  space  bounded  by  Gn  >  a00  (i.e.  axially  sym¬ 
metric  compressive  stresses).  During  the  experiment,  the  strains  en  and  eee  are  meas¬ 
ured  along  with  the  stresses  Gn  and  G00.  Isotropy  of  deformation  in  the  {23}  plane  is 
assumed  (i.e.  e00  =  £22  =  £33)-  Each  of  the  samples  tested  during  the  course  of  this 
work  were  relatively  isotropic,  with  linear  compressibilities  measured  in  different 
orientations  varying  by  less  than  ±5%. 

Intuition  leads  us  to  try  to  formulate  a  rheology  in  terms  of  the  mean  stress  (a), 
the  shear  stress  (x),  the  volumetric  strain  (k),  and  the  shear  strain  (y),  defined  as: 


o  =  (on+  2aee)/3.0 

(la) 

x  =  (gu  -  g00)/2.O 

(lb) 

k  =  £u  +  2e00 

(lc) 

y  =  £n-e00 

(Id) 

For  an  isotropic,  linear  elastic  solid,  K  =  o/K  and  y  =  x/G,  where  K  and  G  are 
the  bulk  and  shear  moduli,  respectively.  Here  we  will  generalize  the  definitions  of  the 
moduli  to  allow  for  the  description  of  nonlinear,  hysteretic  deformation.  We  define  the 
more  general  bulk  and  shear  "moduli"  ( K  and  G)  to  be: 


K  = 

/V 

G  = 


do 

dK 

dr 

dy 


(2a) 

(2b) 


3.4  Pure  Mode  Behavior 

We  begin  by  discussing  the  stress  versus  strain  relationships  for  perturbations  in 
mean  stress  at  constant  shear  stress  and  perturbations  in  shear  stress  at  constant  mean 
stress.  We  refer  to  these  loading  paths  as  "pure"  mode  loading  paths.  Over  a  limited 
but  significant  range  in  mean  stress  perturbations,  it  is  found  that  for  all  rock-types  the 
pure  mode  data  can  be  modeled  with  a  simple  description  of  K  which  will  be  referred 
to  as  the  asymmetric  bow-tie  parameterization  (ABP).  The  model  parameterization,  and 
an  example  comparison  with  data,  is  illustrated  in  Figure  20.  In  this  model,  K  is 
assumed  to  change  linearly  with  mean  stress  for  monotonic  changes  in  mean  stress.  At 
any  reversal  in  the  direction  of  change  of  the  mean  stress,  K  reverts  to  a  value  K', 
which  for  simplicity  here  is  modeled  as  a  linear  function  of  the  mean  stress  at  the 
reversal  in  loading  direction.  It  is  also  observed  that  the  deformation  exhibits  a 
memory  of  the  stress-strain  state  at  previous  extrema  in  mean  stress,  a  characteristic 
commonly  referred  to  as  discrete  memory  (see  Guyer  et  al.  [1995]).  To  enforce  this 
memory  we  require: 


0" 


J 


do 


(3) 
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where  o'  and  o"  are  two  sequential  extrema  in  the  mean  stress  loading  path.  If  we  are 
trying  to  model  k  for  the  loading  path  o'  ->  o"  ->  o'",  we  specify  (or  already  know) 
K}(o)  and  use  equation  (3)  to  find  K2(o).  K2  is  then  used  to  compute  K  during  the 
loading  a"  ->  o'",  and  equation  3  can  then  be  rewritten  to  find  K3  for  the  next  loading 
leg.  With  prior  specification  of  the  functional  form  of  K,  the  necessary  model  parame¬ 
ters  are  constrained  by  the  memory  enforced  by  equation  (3).  For  the  case  of  the 
ABP,  there  are  three  free  parameters,  two  describing  the  mean  stress  dependence  on 
K',  and  one  describing  dK /do  during  the  current  loading  leg. 

For  the  case  of  shear  stress  changes  at  constant  mean  stress,  it  is  found  that  the 
same  parameterization  appears  to  describe  G  (now  in  terms  of  y  and  x).  Thus  the 
ABP  parameterization  appears  to  describe  the  hysteresis  in  both  the  volumetric  and 
shear  strains  (at  constant  shear  and  mean  stress,  respectively).  This  is  found  to  be  true 
for  each  of  the  rock-types  tested. 

Figure  21  illustrates  the  results  of  model  fits  to  data  from  Berea  sandstone  and 
N-tunnel  tuff.  By  comparing  results  from  different  initial  mean  stresses,  it  is  clear  that 
the  model  parameters  describing  K  and  G  are  dependent  on  the  initial  mean  stress 
This  dependence  is  much  stronger  and  functionally  more  complex  for  a  rock  like  Berea 
sandstone  compared  to  tuff.  These  observations  illustrate  that  for  large  perturbations 
in  o,  the  data  require  a  functionally  more  complicated  parameterization,  with  K'  and 
G '  being  described  by  nonlinear  functions  of  o.  Accordingly,  large  amplitude  perturba¬ 
tions  also  require  dK/do  and  dG/dx  during  monotonic  loading  and/or  unloading  to  be 
described  by  higher  order  functions  as  well  (e.g.  see  the  previous  work  on  Young’s 
modulus:  Boitnott  [1993]).  We  should  also  expect  that  K  and  G  may  be  dependent  on 
x  (over  and  above  the  dependency  of  G'  on  x),  and  thus  similar  complexity  should 
arise  for  large  shear  perturbations.  Despite  this  complexity,  over  limited  but  rather 
large  changes  in  o  and  x,  the  ABP  appears  to  provide  a  good  approximation  of  the 
deformation  from  pure  shear  and  pure  mean  stress  perturbations  in  variety  of  rocks  and 
covering  a  wide  range  of  stress  conditions. 

3.5  Some  Notes  on  Notation 

Before  proceeding  with  a  discussion  of  the  mixed  mode  loading  paths,  some  nota- 
tional  details  need  to  be  explained.  Subscripts  are  used  to  distinguish  between  moduli 
measured  from  different  loading  paths.  Loading  paths  and  their  associated  subscripts 
are  as  follows:  mean  stress  (a),  shear  stress  (x),  uniaxial  stress  (us),  and  uniaxial  strain 
(«e).  Elastic  moduli  inferred  from  compressional  and  shear  velocities  (assuming  iso¬ 
tropy  and  linear  elasticity)  are  denoted  by  the  subscript  (v).  Other  common  moduli  are 
also  discussed,  but  unlike  K  and  G ,  are  defined  only  for  specific  loading  paths.  These 
moduli  include  Young’s  modulus  (Eus  =  dan/38n  lAOee=0)  and  the  uniaxial  strain 
modulus  (YuZ  =  5<5i]/d£nlAeee=o)-  1°  similar  fashion,  a  generalized  Poisson’s  ratio 
(9U5  =  -deQe/dE]  1 1 ao60=o)  is  also  used.  Parameters  denoted  without  the  O  represent 
values  inferred  from  other  parameters  assuming  isotropic,  linear  elasticity  (e.g.  Eaz  is 
the  inferred  Young’s  modulus  derived  from  measurements  of  Ka  and  Gx).  The  super¬ 
script  (  ')  is  used  to  denote  the  value  of  a  parameter  measured  immediately  after  a 
reversal  and  loading  direction,  and  symbols  enveloped  in  <  >  denote  the  average  value 


-  12  - 


of  the  enveloped  parameter  over  an  extended  range  of  loads  (such  as  over  an  entire 
loading  and  unloading  cycle). 

3.6  Behavior  During  Mixed  Mode  Loading  Paths 

Unfortunately,  the  data  indicate  that,  in  general,  it  is  not  sufficient  to  superimpose  the 
rheologies  discussed  above  for  Ka  and  Gx  to  describe  deformation  during  arbitrary 
(mixed-mode)  loading  paths.  While  under  some  conditions,  and  for  some  rock  types,  a 
simple  decoupled  model  may  be  appropriate,  additional  effects  are  commonly 
observed.  Figures  22a,  23a  and  24a  plot  the  volumetric  strain  as  a  function  of  the 
mean  stress  for  Berea  sandstone,  NTS  tuff,  and  Sierra  White  granite  during  the  2nd 
perturbation  sequence  (as  shown  in  Figures  18  and  19).  Here  we  see  that  for  the 
mixed  mode  cases  (i.e.  uniaxial  stress  and  uniaxial  strain),  Ka  is  not  sufficient  to 
model  the  volumetric  strain  during  other  modes  of  deformation. 

We  do  however  note  some  pattern  to  the  stress-strain  behavior  shown  in  Figures 
22a,  23a,  and  24a.  In  those  figures  ,  it  can  be  seen  that  the  shear  stress  perturbation 
(at  constant  mean  stress)  causes  an  appreciable  decrease  in  volumetric  strain  (i.e.  a 
dilation).  Accordingly,  the  uniaxial-stress  and  uniaxial-strain  perturbations  produce 
less  volumetric  strain  (per  unit  change  in  mean  stress)  than  is  observed  during  the  pure 
mean  stress  cycles.  This  suggests  that  a  "coupling"  term  is  required  to  handle 
volumetric  strain  resulting  from  changes  in  shear-stress.  Similarly  for  the  shear  strain 
(Figures  22b,  23b,  and  24b),  there  is  again  evidence  of  coupling  in  the  sense  that  the 
shear  stress  versus  shear  strain  relationship  is  not  fully  described  by  Gx  when  a  is 
changing.  Here  again,  a  coupling  term  is  required  to  handle  the  effects  of  mean  stress 
changes  on  the  shear  rheology. 

3.7  Characterizing  the  Coupling  Terms 

The  coupling  terms  appear  to  be  a  strong  function  of  rock  type,  and  for  Berea 
appear  to  become  less  significant  as  the  initial  mean  stress  increases.  It  is  likely  that 
they  result  at  least  in  part  from  stress  induced  anisotropy  in  mechanical  properties  in 
response  to  differential  loads. 

In  many  respects,  the  coupling  terms  appear  relatively  elastic  (i.e.  uni-valued 
functions  of  stress  state),  allowing  for  the  potential  to  develop  a  general  rheology  by 
superposition  of  elastic  coupling  terms  with  hysteretic  and  stress  dependent  Ka  and 
Gv  The  development  of  a  hysteretic  rheology  is  not  fully  completed.  In  its  place,  a 
nonlinear  elastic  (non-hysteretic)  approximation  to  the  data  is  developed  here,  which 
will  prove  useful  as  an  illustration  of  the  general  character  of  the  coupling  terms.  The 
rheology  can  be  described  by: 


+  AtC(o,t) 


(4a) 


Ay  =  -^-  +  Aa-y(a,t) 


(4b) 
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The  terms  £  and  \|/  are  the  coupling  terms.  C,  handles  the  effects  of  changing  shear 
stress  on  the  mean  stress  rheology  and  V|/  handles  the  effects  of  changing  mean  stress 
on  the  shear  rheology. 

While  the  form  of  the  coupling  terms  is  still  not  fully  developed,  reasonable 
behavior  can  be  achieved  using  functionally  similar  forms  for  t,  and  \|/: 


C  =  -Axae~aa 

(4c) 

\| f  =  -Bx^e~ba 

(4d) 

Preliminary  attempts  to  fit  the  Berea  sandstone  data  lead  to  the  following  values: 
<A>  =  1.8e -05  MPa~ia+l),  <a>  ~  1.5,  <a>  ~  0.2  MPa~\  <B>  ~  8.0e-06  MPa~(p+1), 
<(3>  =  1.0,  and  <b>  ~  0.05  MPa'1.  Using  these  values,  along  with  constant  values  of 
<Ka>  =  10  GPa  and  <Gt>  =  8  GPa ,  model  predictions  for  the  data  in  Figure  22  is 
shown  in  Figure  25.  Note  that  the  nonlinear  model  (equations  4a-d)  reproduces  the 
dominant  features  in  the  data,  with  the  exception  of  the  hysteresis  which  in  this  case 
has  been  explicitly  ignored. 

One  underlying  implication  of  the  coupling  terms  is  that  inferred  (average)  elastic 
moduli  are  a  function  of  the  loading  path  and  the  stress  state  from  which  they  are 
derived.  As  an  illustration,  equations  4a  and  4b  can  be  rewritten  as: 


K  = 


Kg 

1  +  rK£ 


(5a) 


G  = 


rGx 


r  +  Gx\)/ 


(5b) 


where  r  is  the  stress  trajectory  parameter  ( r  =  Ax/ Act).  Note  that  if  r  is  positive, 
K  >  KQ  and  G  >  Gx  (since  C,  and  \)/  are  found  to  be  negative). 

3.8  Comparisons  with  dynamic  moduli 

The  coupling  terms  described  by  equations  4c  and  4d  have  been  derived  from 
moderate  stress  amplitude  perturbations,  and  it  remains  unclear  as  to  whether  they 
result  from  stress  hysteretic  phenomenon  or  if  they  are  dominantly  elastic  in  character. 
To  begin  addressing  this  issue  a  second  testing  protocol  was  developed  which  com¬ 
bines  a  complex  loading  history,  static  strain  measurements,  and  ultrasonic  velocity 
measurements.  Compressional  and  shear  velocities  were  measured  using  ultrasonic 
pulse  propagation  (=1  MHz)  along  the  axis  of  the  cores.  As  in  the  previous  experi¬ 
ments,  room  dry  Berea  sandstone  was  used  as  the  sample  material.  The  addition  of 
the  velocity  data  provides  information  concerning  the  mechanical  properties  (i.e. 
moduli)  for  very  small  stress  perturbations,  and  thus  should  provide  constraints  on  the 
behavior  of  C,  and  \|/  in  the  limit  of  small  strains. 
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The  loading  protocol  is  shown  in  Figure  26.  After  an  initial  loading  and  unload¬ 
ing  of  the  confining  pressure  (to  season  the  sample),  confining  pressure  was  increased 
to  50  MPa  and  velocities  were  measured  at  5  points  during  the  loading.  The  confining 
pressure  was  returned  to  5  MPa,  and  a  slight  differential  load  (1  MPa)  was  applied. 
Confining  pressure  was  then  increased  (at  constant  differential  stress),  stopping  at  the 
same  5  points  for  which  velocities  were  measured  during  hydrostatic  loading.  At  each 
of  these  5  points  during  the  loading,  a  series  of  stress  perturbations  was  applied,  con¬ 
sisting  of  perturbations  in  mean  stress,  uniaxial  stress,  and  shear  stress.  Velocities  were 
measured  at  various  points  during  this  loading  sequence. 

The  observed  nonlinearities  in  the  static  stress-strain  relationships  are  quantita¬ 
tively  consistent  with  those  from  the  moderate  strain  experiments  discussed  earlier. 
Stress-strain  relationships  during  the  mean  stress  and  shear  stress  perturbations  are 
used  in  infer  <K0>  and  <GX>.  Stress-strain  relationships  during  the  uniaxial  stress 
perturbations  are  used  in  infer  Young’s  modulus  and  Poisson’s  ratio  (cEus>  and 
<VIW>).  In  addition,  velocities  measured  during  the  hydrostatic  loading  (events  1-5  in 
Figure  26)  are  used,  along  with  the  measured  bulk  density,  to  infer  dynamic  elastic 
constants  controlling  compressional  and  shear  velocities  (Yv  and  Gv).  Assuming  iso¬ 
tropic,  linear  elasticity,  these  three  sets  of  elastic  constants  were  used  to  compute 
(infer)  all  other  elastic  constants.  Thus  there  are  three  sets  of  inferred  elastic  con¬ 
stants,  one  inferred  from  pure  mode  perturbations,  the  second  inferred  from  the  mixed 
mode  perturbation  of  uniaxial  stress,  and  the  third  inferred  from  the  compressional  and 
shear  velocities. 

The  inferred  moduli  are  plotted  in  Figure  27.  Comparing  the  results  of  the  two 
static  measurements,  the  expected  systematic  differences  between  elastic  constants 
inferred  using  pure  modes  and  those  inferred  from  uniaxial  stress  are  clearly  seen. 
Most  notably,  inferred  bulk  modulus  is  higher  for  uniaxial  stress  than  it  is  for  the 
mean  stress  perturbation.  This  is  the  result  of  finite  At Z,  during  the  uniaxial  stress  per¬ 
turbations.  Inferred  shear  modulus  is  not  notably  different  for  the  two  cases,  as  is 
predicted  by  the  nonlinear  rheology  for  these  stress  conditions. 

Comparing  these  results  with  the  moduli  derived  from  the  velocities,  the  statically 
determined  moduli  are  in  general  lower.  Importantly,  the  systematics  between  the  K' s 
and  G’s  with  respect  to  the  mean  stress  are  similar  between  the  pure  mode  static  and 
the  velocity  derived  cases,  and  different  for  the  mixed  mode  static  case  (uniaxial 
stress).  This  difference  in  systematics  with  mean  stress  is  clearly  reflected  in  the  esti¬ 
mates  of  the  v’s,  which  for  isotropic,  linear  elasticity,  is  sensitive  only  to  the  ratio 
K/G. 

3.9  Correcting  for  Strain  Amplitude 

It  has  long  been  realized  that  the  stress  amplitude  of  the  perturbations  used  to 
measure  elastic  constants  strongly  influences  the  measured  values,  and  that  this  is  one 
of  the  dominant  causes  of  discrepancies  between  small  amplitude  dynamic  measure¬ 
ments  and  larger  amplitude  static  measurements  (see  for  example  Walsh.  [1965]). 
Stress  amplitude  effects  arrize  primarily  from  hysteretic  nonlinearities.  In  the  limit  of 
small  strains,  hysteresis  is  minimized  and  it  is  found  that  inferred  moduli  (for  a  partic¬ 
ular  mode  of  deformation)  approach  the  moduli  at  reversals  in  load  observed  during 
larger  amplitude  cycles  (see  Guyer  et  al.  [1997]).  In  order  to  better  compare  the  static 


-  15  - 


and  dynamic  data  collected  here,  the  local  moduli  K'g,  G\,  E'US,  and  V'us  at  each 
reversal  in  loading  direction  were  used  in  place  of  <Ka> ,  <GX>,  <EUS>,  and  <$us>. 
These  stress  amplitude  corrected  moduli  were  then  compared  with  the  dynamic  moduli 
from  the  ultrasonic  velocities.  Results  are  plotted  in  Figure  28.  The  strain  amplitude 
correction  of  the  static  moduli  results  in  an  elevation  of  the  moduli  with  respect  to 
their  predecessors  (i.e.  static  values  in  Figures  27a  and  27b).  Importantly,  the  ampli¬ 
tude  correction  has  not  changed  the  systematics  between  the  K’s,  G’s,  and  v’s  with 
mean  stress  for  either  case.  K'a  and  G\  from  pure  mode  static  measurements,  com¬ 
pare  favorably  with  Kv  and  Gv  from  the  velocities,  although  velocity  derived  moduli 
are  still  slightly  higher.  G'us  is^not  significantly  different  from  G'x,  but  K'us  is  not¬ 
ably  higher  than  the  inferences  K'a. 

3.10  Some  General  Remarks  About  the  Rheology 

The  experiments  clearly  demonstrate  that  the  deformation  of  all  of  the  samples 
tested  is  nonlinear.  The  data,  along  with  a  nonlinear  elastic  model  which  approximates 
much  of  the  observed  nonlinearity,  illustrate  that  elastic  moduli  are  dependent  on  the 
mode  of  deformation  from  which  they  have  been  measured.  The  most  dominant  cou¬ 
pling  effect  in  Berea  sandstone  is  the  C,  term,  which  acts  to  increase  the  bulk  modulus 
(. K )  during  mixed  mode  loading  paths  where  both  mean  stress  and  shear  stress  increase 
(or  decrease)  simultaneously. 

The  magnitude  of  the  coupling  terms  C,  and  \]/  are  found  to  decrease  with  increas¬ 
ing  mean  stress,  and  increase  with  increasing  shear  stress.  They  likely  reflect,  at  least 
in  part,  shear  stress  (i.e.  differential  stress)  induced  anisotropy.  As  differential  stress  is 
increased,  the  stiffness  to  axial  deformation  is  increased,  thus  resisting  both  volumetric 
and  shear  strains.  In  addition,  the  coupling  terms  may  also  reflect  shear  stress  induced 
dilatancy,  as  grain  and  crack  surfaces  attempt  to  slide  past  opposing  features. 

There  is  some  limited  evidence  that  the  stress  induced  anisotropy  mechanism  is 
the  dominant  process,  as  the  coupling  terms  controlling  moderate  stress  amplitude 
moduli  (e.g.  <K>  and  <G>)  also  appear  to  control  strain  amplitude  corrected  moduli 
(e.g.  K'  and  G').  Evidence  of  shear  induced  dilatancy  is  less  clear,  but  is  appealing  in 
that  dilatancy  is  known  to  occur  at  higher  differential  stresses.  It  is  interesting  to  note 
that  the  nonlinear  rheology  presented  in  equations  4a-d  predicts  dilatancy  (i.e.  A'cO.O) 
at  high  differential  stresses.  While  extension  to  higher  differential  stresses  will  require 
higher  order  functions  for  £  and  \\f  (as  well  as  hysteresis  in  these  terms),  it  is  reason¬ 
able  to  expect  that  the  rheology  is  extendible  to  these  conditions. 

Development  of  a  hysteretic  rheology  which  includes  all  of  the  observed  effects 
of  loading  path  still  requires  more  work.  Reasonable  behavior  can  be  attained  by  using 
hysteretic  pure  mode  moduli  ( Ka  and  GT)  and  combining  them  with  nonhysteretic  cou¬ 
pling  terms  (i.e.  placing  hysteretic  Ka  and  Gx  in  equations  4a, b).  The  model  does  not 
produce  the  observed  reduction  in  hysteresis  of  K  during  uniaxial  stress  cycles  and 
also  tends  to  over-predict  hysteresis  in  G  during  mixed  mode  loading.  In  the  context 
of  the  rheology  developed  here,  these  effects  require  either  the  introduction  of  hys¬ 
teretic  coupling  terms  or  the  replacement  of  Ka  and  Gx  in  equations  4a, b  with  loading 
path  dependent  terms.  It  should  also  be  noted  that  the  coupling  terms  described  here 
could  easily  be  introduced  into  the  endochronic  rheology  of  Xu  et  al.  [1996].  Their 
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endochronic  rheology  is  essentially  nothing  more  than  a  simple  decoupled  superposi¬ 
tion  of  Kg  and  GT  (using  a  different  parameterization  than  developed  here),  and  thus 
adding  elastic  coupling  terms  should  be  straightforward.  Despite  the  shortcomings  of 
the  use  of  nonhysteretic  coupling  terms,  the  resulting  rheology  mimics  observation 
remarkably  well,  and  thus  should  be  suitable  to  assess  whether  the  observed  nonlinear¬ 
ities  are  important  in  the  context  of  a  particular  application. 

4.  THE  EXPLOSION  LOADING  PATH 

4.1  A  Simpler  Rheology  for  an  Important  Loading  Path 

Despite  the  apparent  complexity  introduced  by  the  coupling  terms,  we  find  that  a 
related  but  simpler  rheology  can  be  formulated  which  is  applicable  to  the  case  of 
compressional  wave  propagation  from  an  explosive  source.  Figure  29  illustrates  the 
loading  protocol  of  an  experiment  designed  to  characterize  deformation  during  loading 
paths  which  mimic  (quasi-statically)  the  stress  path  expected  as  a  result  of  an  under¬ 
ground  explosion.  The  loading  paths  are  constructed  to  follow  that  predicted  for  a 
Heaviside  step  function  in  radial  stress  at  a  spherical  cavity  in  a  linear  elastic  material. 
The  resulting  explosion  loading  path  involves  a  uniaxial  strain  loading  followed  by  a 
mixed-mode  unloading  which  follows  a  damped  oscillation  about  a  final  stress  state. 
The  unloading  path  is  a  function  of  the  range  (distance  from  source)  and  the  effective 
elastic  constants  of  the  material.  In  the  experiment  illustrated  in  Figure  29,  we  have 
performed  four  cycles  of  uniaxial  strain  loading,  followed  by  three  successive  explo¬ 
sion  loading  paths  with  increasing  amplitude  (representing  decreasing  range). 

Figure  30  illustrates  the  stress-strain  relationship  during  this  experiment,  decon¬ 
volved  into  stresses  and  strains  as  defined  in  equations  l[a,b,c,d].  Here  we  see  that  the 
deformation  is  well  behaved  (although  nonlinear  and  hysteretic),  without  noticeable 
complications  by  the  coupling  terms  (as  was  seen  in  Figure  22).  In  fact,  the  data  can 
be  modeled  fairly  well  using  a  uncoupled  ABP  for  both  Kue  and  Gue,  (i.e.  where  K 
and  G  are  constrained  by  the  uniaxial  strain  data).  Note  how  the  loops  for  the  explo¬ 
sion  paths  are  not  easily  distinguishable  from  the  uniaxial  strain  loops  which  preceded 
them.  In  effect,  this  short-cut  groups  the  coupling  terms  discussed  previousing  into  the 
ABP  description  of  Kue  and  GuV  and  is  thus  valid  only  for  loading  paths  for  which  the 
coupling  terms  themselves  are  functionally  consistent  with  the  ABP.  The  explosion 
loading  path  examined  here  appears  to  be  such  a  path,  and  thus  this  simplified  model 
appears  appropriate  for  the  problem  of  the  spherically  propagating  compressional  wave 
for  an  explosion-like  source. 

Figure  31  illustrates  the  stress-strain  relationship  for  N-Tunnel  tuff  during  an 
experiment  of  the  type  illustrated  in  Figure  29.  Again,  the  deformation  (both 
volumetric  and  shear)  for  the  explosion  loading  path  is  found  to  be  modelable  using  an 
"uncoupled"  rheology  constrained  by  the  uniaxial  strain  path. 

In  Figure  32,  we  compare  the  measured  strain  paths  with  the  strain  paths 
predicted  by  a  linear  elastic  rheology  (for  the  identical  stress  paths)  for  two  explosion 
paths  in  Berea  Sandstone.  The  elastic  constants  used  for  the  linear  elastic  predictions 
are  constrained  to  enforce  the  uniaxial  strain  loading  leg  during  the  early  stages  of  the 
loading  path.  The  measured  strain  path  upon  unloading  is  noticeably  different  from 
that  of  the  linear  elastic  model.  Such  differences  are  readily  noted  even  for  conditions 
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which  would  typically  be  considered  well  outside  the  elastic  radius. 

The  implications  of  these  differences  are  not  well  known  at  this  point.  However, 
qualitatively  one  can  reason  that  if  the  measured  strains  were  used  to  infer  a  stress 
path  (using  a  linear  elastic  model),  the  inferred  stress  path  would  be  different  from  the 
actual  stress  path.  Thus,  if  back-propagated  to  the  source,  the  error  would  result  in  an 
incorrect  (unphysical)  source  function.  In  addition  to  this  contamination  of  an  observa¬ 
tion  at  a  particular  range,  the  rheology  will  also  have  effects  on  the  propagation  of  the 
compressional  wave  itself,  further  distorting  a  source  function  derived  through  inap¬ 
propriate  use  of  a  linear  elastic  model. 

4.2  Observations  at  Larger  Stress  Amplitudes 

A  few  exploratory  tests  have  been  performed  in  an  attempt  to  see  when  the  char¬ 
acter  of  the  deformation  changes  as  we  go  to  higher  stress  amplitudes.  At  an  initial 
mean  stress  of  10  MPa,  an  experiment  on  Berea  sandstone  (Figure  33)  with  peak  axial 
stresses  of  80  MPa  illustrates  that  the  same  rheology  appears  to  hold  even  for  these 
large  strains.  A  experiment  on  N-Tunnel  tuff  behaved  similarly  for  peak  axial  stresses 
of  40  MPa,  but  caused  failure  during  the  uniaxial  strain  loading  during  an  attempt  to 
go  to  80  MPa.  These  experiments  indicate  that  the  rheology  appears  viable  over  the 
entire  moderate  strain  regime  (until  the  onset  of  damage). 

In  the  small  strain  limit,  the  rheology  approaches  a  non-linear  elastic  model  at 
low  mean  stresses  and  a  linear  elastic  model  at  high  ambient  mean  stresses.  In  previ¬ 
ous  work,  the  rheology  has  been  shown  to  remain  viable  to  strain  amplitudes  as  low  a 
10~5. 


5.  CONFINED  COMPRESSION  AND  DAMAGE 


Experiments  have  also  been  conducted  on  Barre  granite  as  exploratory  experi¬ 
ments  for  studying  the  effect  of  confining  pressure  on  uniaxial  stress  deformation  and 
for  extending  the  moderate  strain  regime  rheology  into  the  brittle  damage  regime.  Fig¬ 
ure  34  illustrates  and  example  of  the  stress  history  of  one  of  the  the  experiments. 
Experiments  have  been  performed  at  two  different  confining  pressures  (0  MPa  and  10 
MPa  respectively). 

Stress-strain  results  a  shown  in  Figure  35a.  During  the  loading,  numerous  small 
perturbations  in  the  load  were  imposed  to  measure  the  unloading  stiffness.  In  these 
experiments  we  find  that  once  we  enter  the  damage  regime  (i.e.  differential  stress  lev¬ 
els  above  150  MPa),  the  stress-strain  relationship  during  cyclic  loading  exhibits  a 
ratcheting  effect,  which  becomes  more  pronounced  as  we  approach  failure.  This 
ratcheting  is  thought  to  reflect  time  dependent  crack  growth,  but  may  also  include 
static  (rate  independent)  phenomena  as  well. 

The  data  provides  some  preliminary  constraint  on  confining  pressure  effects  in  the 
moderate-strain  and  damage  regimes.  In  Figure  35b,  The  logarithm  of  the  local  Youngs 
modulus  ( E )  during  loading  for  each  experiment  is  plotted  using  solid  lines  and  the 
logarithm  of  the  unloading  stiffnesses  ( E ')  are  plotted  using  circles.  Looking  at  the 
data  in  the  moderate-strain  regime  (differential  stresses  below  150  MPa),  we  see  that 
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the  effect  of  confining  pressure  on  both  E  and  E'  can  be  modeled  with  a  simple  scale 
factor  which  is  of  comparable  magnitude  for  both  E  and  E' .  This  is  consistent  with 
the  X  parameter  as  discussed  earlier.  As  discussed  in  the  context  of  the  uniaxial  stress 
data  in  section  2,  the  data  here  indicates  that  X  is  a  function  of  differential  stress  at 
low  axial  loads.  Also  note  that  once  in  the  brittle-damage  regime,  this  simple  scaling 
with  confining  pressure  breaks  down,  as  the  effect  of  confining  pressure  on  E  develops 
a  more  pronounced  pressure  effect  than  on  E' ,  which  is  likely  due  in  large  part  to  time 
dependent  damage  and  resulting  deformation. 
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Radial  Strain  (xl.Oe+06) 


Figure  1:  Axial  stress  versus  axial  strain  and  axial  stress  versus  radial  strain  for  cyclic 
loading  of  Sierra  White  granite  in  unconfined  compression.  Note  the  absence  of 
hysteresis  in  the  radial  strain.  This  observation  allows  for  the  possibility  of  a  sim¬ 
ple  model  incorporating  the  effects  of  confining  pressure  on  the  hysteresis  in  axial 
strain. 
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Figure  2:  Comparison  of  Experimental  data  (from  Hilbert  et  al.,  1994)  with  model 
predictions  for  hysteretic  deformation  in  Berea  Sandstone.  A  simple  scaling  law  is 
used  to  incorporate  the  effects  of  changing  confining  pressure  on  the  axial 
stiffness  during  uniaxial  strain  loading.  The  simple  model  appears  to  contain 
many  of  the  features  exhibited  by  the  data.  Model  results  are  based  on  experi¬ 
ments  using  a  different  sample  of  Berea  sandstone,  and  thus  small  quantitative 
discrepancies  between  model  and  data  should  not  to  be  considered  significant. 
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Figure  3:  Example  of  the  loading  protocol  used  to  constrain  the  effect  of  confining 
pressure  on  the  uniaxial  stress  rheology.  At  a  series  of  different  confining  pres¬ 
sures  (Pc  =  5,  10,  and  20  MPa)  A  sequence  of  axial  stress  perturbations  were 
made  while  axial  and  radial  strain  were  monitored. 


Figure  4:  Axial  stress  versus  axial  strain  for  Berea  sandstone  during  the  axial  loading 
sequences  at  5  and  20  MPa  confining  pressure  shown  in  figure  3. 


0 


Axial  Stress  (MPa) 


20 


Figure  5:  Local  Young’s  modulus  (E  =  3an/9en)  for  Berea  sandstone  during  uniaxial 
stress  cycles  at  5  and  20  MPa  confining  pressure.  Data  is  from  the  second  large 
amplitude  stress  cycle  at  each  pressure  in  the  sequence  shown  in  Figure  3. 
Arrows  indicate  loading  direction. 
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Figure  6:  Local  Young’s  modulus  ( E  =3a11/3e11)  for  N-Tunnel  tuff  during  uniaxial 
stress  cycles  at  5  and  20  MPa  confining  pressure.  Data  is  from  the  second  large 
amplitude  stress  cycle  at  each  pressure  in  the  sequence  similar  to  that  shown  in 
Figure  3.  Arrows  indicate  loading  direction. 
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Sierra  White  Granite 


Figure  7;  Local  Young’s  modulus  ( E  —  3<Tjj/3£h)  for  Sierra  White  granite  during 
uniaxial  stress  cycles  at  5  and  20  MPa  confining  pressure.  Data  is  from  the 
second  large  amplitude  stress  cycle  at  each  pressure  in  the  sequence  similar  to 
that  shown  in  Figure  3.  Arrows  indicate  loading  direction. 
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Figure  8:  Stress  versus  strain  for  Berea  sandstone  during  a  uniaxial  stress  cycle  at  a 
confining  pressure  of  10  MPa.  a)  Axial  stress  versus  axial  strain,  b)  Radial  stress 
versus  radial  strain,  c)  Mean  stress  versus  volumetric  strain,  d)  Shear  stress 
versus  shear  strain.  Note  the  strong  hysteresis  in  (b)  and  the  lack  of  hysteresis  in 
(c).  In  this  case  the  hysteresis  can  be  thought  to  be  confined  to  the  shear  rheol- 
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Figure  9:  Stress  versus  strain  for  N-Tunnel  tuff  during  a  uniaxial  stress  cycle  at  a 
confining  pressure  of  10  MPa.  a)  Axial  stress  versus  axial  strain,  b)  Radial  stress 
versus  radial  strain,  c)  Mean  stress  versus  volumetric  strain,  d)  Shear  stress 
versus  shear  strain.  Note  the  measurable  hysteresis  in  each  plot,  indicating  that 
both  the  mean  stress  and  the  shear  stress  rheologies  are  hysteretic. 
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Figure  10:  Stress  versus  strain  for  Sierra  White  granite  during  a  uniaxial  stress  cycle  at 
a  confining  pressure  of  10  MPa.  a)  Axial  stress  versus  axial  strain,  b)  Radial 
stress  versus  radial  strain,  c)  Mean  stress  versus  volumetric  strain,  d)  Shear  stress 
versus  shear  strain.  Note  the  measurable  hysteresis  in  (b)  which  was  not 
observed  during  the  unconfined  test  in  figure  1 . 
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Figure  12:  Volumetric  strain  versus  confining  pressure  (mean  stress)  for  Berea  sand¬ 
stone  durng  the  loading  sequence  in  figure  11.  Note  the  strong  hysteresis  here  in 
contrast  to  the  lack  of  hysteresis  in  the  mean  stress  versus  volumetric  strain  rela¬ 
tionship  during  uniaxial  stress  cycling  (compare  with  figure  8c). 
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Figure  14:  Volumetric  strain  versus  confining  pressure  (mean  stress)  for  Sierra  White 
granite  durng  the  loading  sequence  in  figure  1 1 .  Note  the  strong  hysteresis  here  in 
contrast  to  the  lack  of  hysteresis  in  the  mean  stress  versus  volumetric  strain  rela¬ 
tionship  during  uniaxial  stress  cycling  (compare  with  figure  10c). 
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Figure  15:  Measured  shear  stress  versus  shear  strain  for  Berea  sandstone.  In  the  left 
figure,  the  effect  of  shear  stress  oscillation  amplitude  is  shown  at  constant  axial 
load  (5.4  MPa).  Note  that  the  smaller  loop  is  nested  within  and  rotated  with 
respect  to  the  larger  loop.  In  the  right  figure,  the  effect  of  normal  load  is  illus¬ 
trated,  with  loops  at  three  different  normal  loads  (3.6,  5.4,  and  7.2  MPa).  Note 
that  the  shear  modulus  increase  and  the  loop  width  decreases  with  increasing  nor¬ 
mal  load. 
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Figure  16:  An  illustration  of  the  simple  "bow-tie"  rheology  used  to  model  the  data.  In 
the  "bow-tie"  example,  the  shear  modulus  is  assumed  to  follow  the  form 
G  =  G'  -  A-\x-xr  I,  where  G’  is  a  function  of  the  axial  load  (On),  x  is  the  shear 
stress,  and  xr  is  the  shear  stress  at  the  last  reversal  in  loading  direction. 
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Figure  18:  Loading  protocol  of  an  experiment  designed  to  characterize  the  nonlinear, 
hysteretic  rheology  of  rock.  The  protocol  involves  repeated  execution  of  a 
sequence  of  loading  modes,  including  perturbations  in  mean  stress,  uniaxial  stress, 
shear  stress,  and  uniaxial  strain.  See  detail  in  Figure  19. 
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Figure  19:  Detail  of  the  loading  sequence  used  in  the  experiment  in  Figure  18.  Ampli¬ 
tudes  of  the  various  perturbations  produce  similar  deviations  in  the  shear  and 
mean  stresses. 
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Figure  20:  Comparison  of  stress-strain  data  and  model  predictions  for  a  mean  stress 
perturbation  on  Berea  sandstone  (left  figure).  The  ABP  parameterization  is  shown 
graphically  in  the  right  figure.  The  local  modulus  ( K )  changes  linearly  with  stress 
during  monotonic  loading  and/or  unloading.  At  each  reversal  in  loading  direction, 
the  local  modulus  reverts  to  a  value  of  K',  which  is  again  a  linear  function  of 
stress.  The  model  is  constrained  by  equation  3  to  enforce  discrete  memory. 
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Figure  21:  Model  results  for  Berea  sandstone  and  NTS  tuff  illustrating  the  dependence 
of  model  parameters  on  mean  stress.  Models  results  for  mean  stress  perturbations 
are  shown  in  the  left  figure  (for  a  detailed  example  see  Figure  20).  Fits  are  to 
data  from  the  mean  stress  perturbations  during  the  1st  and  3rd  sequence  in  Figure 
18.  The  other  figures  (center  ->  Berea;  right  ->  Tuff)  illustrate  the  inferred  shear 
rheology,  again  for  the  1st  and  3rd  sequence  in  Figure  18.  While  the  ABP  holds 
for  moderate  perturbations  in  stress,  large  perturbations  require  a  functionally 
more  complex  parameterization.  This  is  particularly  true  for  Berea  sandstone, 
which  exhibits  a  large  mean  stress  dependence. 
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Figure  22:  Volumetric  strain  versus  mean  stress  and  shear  strain  versus  shear  stress  for 
Berea  during  the  second  loading  sequence  in  Figure  18.  Note  that  the  shear  stress 
produces  a  significant  dilation.  Note  also  that  the  pure  mode  stress  strain 
response  is  not  reproduced  during  mixed  mode  loading  cycles  for  either  the  shear 
or  the  mean  stress  cases. 
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Figure  23:  Volumetric  strain  versus  mean  stress  and  shear  strain  versus  shear  stress  for 
N-Tunnel  tuff  during  the  second  loading  sequence  in  Figure  18.  Note  that  the 
stress  versus  strain  response  is  qualitatively  similar  to  that  of  Berea  (compare 
with  Figure  22),  but  that  many  of  the  nonlinear  effects  are  not  as  exaggerated. 
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Figure  24:  Volumetric  strain  versus  mean  stress  and  shear  strain  versus  shear  stress  for 
Sierra  White  granite  during  the  second  loading  sequence  in  Figure  18.  Note  that 
the  stress  versus  strain  response  is  qualitatively  similar  to  that  of  Berea  (compare 
with  Figure  22),  but  that  many  of  the  nonlinear  effects  are  not  as  exaggerated. 


Volumetric  Strain 


-  41  - 


Figure  25:  Predictions  of  the  nonlinear  elastic  model  for  the  data  in  Figure  22.  Note 
the  general  agreement  between  the  model  (equations  4a-d)  and  the  data  for  the 
entire  loading  sequence. 
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Figure  26:  Loading  protocol  for  an  experiment  designed  to  integrate  the  moderate 
strain  rheology  with  ultrasonic  velocity  data.  Following  two  cycles  in  hydrostatic 
confining  stress,  as  sequence  of  perturbations  in  stress  is  repeated  five  times  from 
five  different  starting  stresses.  Perturbation  sequences  consisted  of  cycles  in  mean 
stress,  uniaxial  stress,  and  shear  stress.  Numbers  indicates  points  in  the  loading 
sequence  where  ultrasonic  velocities  were  measured. 
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Figure  27:  Inferred  ^elastic  constants  determined  from:  a)  [top  left]  pure  mode 
measurements  of  <Ka>  and  <GT>;  b)  [top  right]  measurements  of  <EUS>  and 
<QUS>  from  uniaxial  stress  perturbations;  and  c)  [bottom  center]  ultrasonic 
compressional  and  shear  wave  velocities  under  hydrostatic  compression  (i.e. 
points  1-5  in  Figure  26). 
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Figure  28:  Inferred  elastic  constants  determined  from:  (a)  [top^pure  mode  meas- 
urements  of  K'a  and  G\;  and  (b)  [bottom]  measurements  of  E'us  and  Vus  from 
uniaxial  stress  perturbations. 
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Figure  29:  Loading  protocol  during  an  experiment  designed  to  characterize  a  simplified 
rheology  for  the  case  of  a  spherically  propagating  compressional  wave.  Uniaxial 
strain  perturbations  are  followed  by  three  "explosion  like"  loading  paths  of 
increasing  amplitude. 
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Figure  30:  Volumetric  strain  versus  mean  stress  and  shear  strain  versus  shear  stress  for 
Berea  sandstone  during  the  loading  protocol  in  Figure  29.  Note  the  well  behaved 
nature  of  the  hysteresis  and  nonlinearity,  exhibiting  the  characteristic  shape  and 
"nested  loop"  behavior  of  the  ABP. 
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Figure  31  Volumetric  strain  versus  mean  stress  and  shear  strain  versus  shear  stress  for 
N-Tunnel  tuff  during  the  loading  protocol  in  Figure  29.  Note  the  well  behaved 
nature  of  the  hysteresis  and  nonlinearity,  exhibiting  the  characteristic  shape  and 
"nested  loop"  behavior  of  the  ABP. 
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Figure  32:  Strain  path  measured  during  the  2nd  and  3rd  "explosions"  as  shown  in  Fig¬ 
ure  6.  The  predictions  for  a  linear  elastic  solid  (with  elastic  constants  constrained 
to  fit  the  uniaxial  strain  loading),  are  shown  for  comparison.  While  qualitatively 
similar,  systematic  differences  arise  from  the  nonlinearity  and  hysteresis  which 
may  have  a  significant  effect  on  wave  propagation  and/or  the  interpretation  of  a 
source  function. 
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Figure  33:  Stress  versus  strain  relationship  for  Berea  sandstone  during  uniaxial  strain 
cycling  similar  to  that  in  Figure  29,  but  with  peak  stress  amplitudes  4  times 
higher.  Note  the  similarity  with  Figure  30,  indicating  that  the  rheology  does  not 
change  even  at  considerably  higher  strain  amplitude  conditions. 
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Figure  34:  Loading  protocol  for  an  experiment  designed  to  probe  the  transition  into  the 
damage  regime.  At  a  series  of  increasing  levels  of  axial  load  (at  fixed  confining 
pressure),  small  perturbations  in  axial  load  are  performed  to  characterize  the 
unloading  modulus  and  its  time  (and  cycle)  dependence.  The  axial  load  is 
increased  until  failure  occurs. 
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Figure  35:  a)  Stress  strain  data  from  two  uniaxial  compression  experiments  on  Barre 
granite,  one  unconfined  and  the  other  at  a  fixed  confining  pressure  of  10  MPa. 
During  the  loading,  perturbations  in  axial  load  were  performed  to  measure  the 
unloading  stiffness,  b)  The  logarithm  of  the  Youngs  modulus  plotted  as  a  function 
of  differential  stress.  Solid  lines  indicating  the  modulus  upon  loading  and  the  cir¬ 
cles  indicate  the  average  modulus  during  the  small  cyclic  perturbations. 


